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ABSTRACT 

Context. Mass loss is essential for massive star evolution, thus also for the variety of astrophysical applications relying on its predic- 
tions. However, mass-loss rates currently in use for hot, massive stars have recently been seriously questioned, mainly because of the 
effects of wind clumping. 

Aims. We investigate the impact of clumping on diagnostic ultraviolet resonance and optical recombination lines often used to derive 
empirical mass-loss rates of hot stars. Optically thick clumps, a non-void interclump medium, and a non-monotonic velocity field are 
all accounted for in a single model. The line formation is first theoretically studied, after which an exemplary multi-diagnostic study 
of an O-supergiant is performed. 

Methods. We used 2D and 3D stochastic and radiation-hydrodynamic wind models, constructed by assembling ID snapshots in radi- 
ally independent slices. To compute synthetic spectra, we developed and used detailed radiative transfer codes for both recombination 
lines (solving the 'formal integral') and resonance lines (using a Monte-Carlo approach). In addition, we propose an analytic method 
to model these lines in clumpy winds, which does not rely on optically thin clumping. 

Results. The importance of the 'vorosity' effect for line formation in clumpy winds is emphasized. Resonance lines are generally 
more affected by optically thick clumping than recombination lines. Synthetic spectra calculated directly from current radiation- 
hydrodynamic wind models of the line-driven instability are unable to in parallel reproduce strategic optical and ultraviolet lines for 
the Galactic O-supergiant A Cep. Using our stochastic wind models, we obtain consistent fits essentially by increasing the clumping 
in the inner wind. A mass-loss rate is derived that is approximately two times lower than what is predicted by the line-driven wind 
theory, but much higher than the corresponding rate derived when assuming optically thin clumps. Our analytic formulation for line 
formation is used to demonstrate the potential importance of optically thick clumping in diagnostic lines in so-called weak-winded 
stars and to confirm recent results that resonance doublets may be used as tracers of wind structure and optically thick clumping. 
Conclusions. We confirm earlier results that a re-investigation of the structures in the inner wind predicted by line-driven instability 
simulations is needed. Our derived mass-loss rate for A Cep suggests that only moderate reductions of current mass-loss predictions 
for OB-stars are necessary, but this nevertheless prompts investigations on feedback effects from optically thick clumping on the 
steady-state, NLTE wind models used for quantitative spectroscopy. 

Key words, stars: early-type - stars: mass-loss - radiative transfer - line: formation - hydrodynamics - instabilities 

1. Introduction theoretical as well as observational evidence of an inhomoge- 

neous, time-dependent wind has become ov erwhelming in th e 

.... e A . ,• c e A . past years (for a comprehensive summary, see Puis et al. 2008). 
Massive stars are fundamental in many fields of modern astro- F 1 F 1 al 

physics. In the present Universe, they dynamically and chem- 
ically shape their surroundings and the interstellar medium by Direct simulations of the time-dependent wind have con- 
their output of ionizing radiation, energy and momentum, and firmed that the so-called line-driven instability c auses a highly 
nuclear processed material. In the distant Universe, they domi- structured wind in bo t h density and velocity (lOwocki et alJ 
nate the ultraviolet (UV) light from young galaxies. Indeed, mas- 1988; iFeldmeierf [19951: iDessart & Owockil [2005). Much indi- 
sive s tars may be regarded as 'cosmic engines' (Bresolin et afj rect evidence of such small-scale inhomogeneities (clumping) 
120081) . Hot, massive stars possess strong and powerful winds that has arisen from quantitative spectroscopy. Clumping has severe 
affect evolutionary time scales, chemical surface abundances, consequences for any interpretation of observed spectra, with 
and luminosities. In fact, changing the mass-loss rates of mas- the inferred mass-loss rates particularly affected. When deriving 
sive stars by only a factor of two has a dramatic effect on their mass-loss rates from observations, wind clumping has tradition- 
overall evolution ( Mev net et al J 1 1 9941) . The winds from these ally been accounted for by assuming optically thin clumps and a 
stars are described by the radiative line-driven wind theory, void interclump medium, while keeping a smooth velocity field, 
in which the standard model (based on the pioneering works Results based on this microclumping approach have, for exam- 
by iLucv & Solomonlll970t ICastor et al.lll975|) assumes that the pie, led to a downward revision of empirical mass-loss rates from 
wind is stationary, spherically symmetric, a nd homogeneous . Wolf-Rayet (W R) stars by roughly a factor of three (reviewed in 
Despite this theory's apparent success (e.g.. IVinketaLl l2000h . ICrowtheill2007h . 
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However, for O stars, highly clumped winds with very low 
mass-loss rates must be invoked in order to reconcile investiga- 
tions of different diagnostics within the microclumping model. 
The most a larming examp l e was the phosphorus v (Pv) UV 
analysis by Fullerton et al. (2006), which indicated reductions 
of previously accepted values by an order of magnitude (or even 
more ), with dwarfs, giants , and supergiants all affected (but see 
also IWaldron & Cassinellil 2010, who argued that xuv radiation 
could seriously alter the ionization fractions of Pv). Such low 
mass-loss rates would be in stark contrast with the predictions 
of line-driven wind theory and have dramatic consequences for 
the evolution of, and feedback from, massive stars. Naturally, 
the widely discrepant values inferred from different observa- 
tions and diagnostics drastically lower the reliability of mass- 
loss rates currently in use, and an explanation is urgently needed. 
A key question is whether the microclumping model fails to de- 
liver accurate empirical rates under certain conditions. 



Simplified techniques to account for opticall y thick clumps 
in X-ray line formation h ave been developed dFeldmeier et al.l 
120031: lOwocki et aT1l2004l) . but it has yet to be settled whether 
or not this is important to consider wh en deriving empirical 
mass-loss rates fr om these diagnostics (Oskinova et al. 2006; 
ICohen et al.l 120 10). First attempts to relax the assumptions of 
th e microclumping mode l for UV resonance line s were made 
bvlOskinova et al.l £2007) (optically thick clumps), Zsarg o et al.l 
(2008) (a non-void interclump medium), a nd Owocki (2008) 
(a non-monotonic velocity field). Sundqvis t et al.l d2010l) (here- 
after Paper I) carried out the first detailed investigation, relaxing 
all the above assumptions, and showed that, indeed, the micro- 
clumping approximation is not a suitable assumption for UV 
resonance line formation under conditions prevailing in typi- 
cal OB-star winds. Recently, these resul ts were empirica l ly sup - 
ported for the case of B supergiants by Prinia & Massa (2010), 
who analyzed profile-strength ratios of the individual compo- 
nents of resonance line doublets and found that the observed 
ratios were inconsistent with lines formed in a smooth or 'mi- 
croclumped' wind. Furthermore, Paper I demonstrated that res- 
onance line profiles calculated from 2D, stochastic wind mod- 
els were compatible with mass-loss rates an order of magnitude 
higher than those derived from the same lines but using the mi- 
croclumping technique. However, as pointed out in that paper, 
a consistent modeling of the resonance lines also introduces de- 
generacies among the parameters used to define the wind struc- 
ture, degeneracies that can only be broken by considering differ- 
ent diagnostics (depending on different parameters) in parallel. 



Here we make a first attempt toward such multi-diagnostic 
studies. We extend our 2D wind models from PapeiT to 3D, 
and relax the microclumping approximation also for the optical 
mass loss diagnostics H a and He n 4686 A (Sect. [2j. In Sect. [3] 
we theoretically investigate H a and resonance line formation in 
clumpy winds, and propose an analytic treatment of the lines that 
does not rely on the microclumping approximation. A simulta- 
neous optical and UV diagnostic analysis is carried out in Sect. [4] 
for the Galactic 06 supergiant A Cep, using time-dependent 
radiation-hydrodynamic (RH) models as well as stochastic ones 
together with our new tools for the radiative transfer in clumped 
winds. These results are discussed in Sect. [5] while two initial 
applications of our analytic formulation are given in Sect. [6] We 
summarize the paper and outline future work in Sect. [7] 



Table 1. Parameters for the time-dependent RH model of 
A Cep (see text). 



Name 


Parameter 


Value 


Spectral type 




06 I(n)fp 


Effective temperature 


Teg 


36 000 K 


Stellar radius 




21.1 R Q 


Surface gravity 


logg 


3.55 


Stellar luminosity 


log L/L 


5.83 


Terminal speed 


Ceo 


2200 km s" 1 


Mass-loss rate 


M 


1.5 x lO^Mo/yr 


Helium abundance 


*He = «Hc/«H 


0.1 


CAK exponent 


Q"CAK 


0.7 


Initial Langevin 


Cturb / Csound 


0.5 


turbulence fluctuation 







Table 2. Basic structure parameters defining a stochastic wind 
model. 



Name 


Parameter 


Clumping factor" 


/c! 


Average time interval 


St 


between release of clumps 




Interclump medium density parameter 


-Me 


Velocity span of clump 


Sv 



" fd may be replaced by the volume filling factor / v . 
The two are related via x ic (see Paper I). 



2. Wind models and radiative transfer 

We create 2D and 3D RH and stochastic wind models by as- 
sembling snapshots in radially independent wind slices (see 
SectO. 

Table 1 summarizes properties of a time -dependent RH 
model computed following the approach in Feldmeier et al.l 
dl997h . which introduces base perturbations from Langevin tur- 
bulence into an unstable line-driven wind. The line force is 
computed with the n onlocal 'Smooth Source Function' (SSF; 
Owocki & Puil l996) method that allows one to follow the non- 
linear evolution of the strong, intrinsic line-deshadowi ng in- 
stabil ity, while also accounting for the diffuse line drag dLucvl 
1984) that reduces (and even eliminates) the instability near the 
wind base. The net result is a highly structured wind character- 
ized by high-speed rarefactions and slower, dense clumps (ac- 
tually shells in these 1-D simulati ons). In comparison to self - 
excited instability simulations (e.g jRunacres & O wocki 2002), 
the base perturbations here induce a somewhat lower onset and 
greater velocity dispersion of the wind clumping. A central goal 
here is to examine the effects of this extensive structure on 
wind diagnost i cs. Ste llar and wind parameters are taken from 
Reno lust et al.l (120041) . except for the mass-loss rate (see Sect. 
4). 

Basic assumptions of our empirical, stochastic models were 
described in detail in Paper I. Essentially, they are constructed 
so to resemble the main structures predicted by the RH simula- 
tions, while still allowing for a variation in the key parameters 
controlling the line formation (see below). 

2. 1 . Parameters describing a structured wind 

When creating our stochastic wind models, we take an heuris- 
tic approach and use a set of parameters to define the structured 
medium. The clumping factor f c \(v)= (p 2 )/(p) 2 , with the an- 
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gle brackets denoting spatially averaged quantities, is the only 
necessary structure parameter when calculating spectra via the 
microclumping technique. Microclumping gives rise to the well 
known result that the opacities for processes that depend on the 
square of the density (for example H^ emission in OB-stars) are 
augmented by / c i as compared to a smooth model with the same 
mass-loss rate; in contrast, opacities for processes that depend 
linearly on the density (for example the uv resonance lines) are 
not directly affected. Thus, if the wind is clumped, mass-loss 
rates derived from smooth models applied to H a are overesti- 
mated by a factor of Jfa. In addition, the occupation numbers 
are modified for all diagnostics because of the changed rates in 
the statistical equilibrium equations. For a comprehensive dis- 
cus siqnon_fhe effects of microclumping on various diagnostics, 
see iPuls etal](l2008l) . 

If the assumptions behind the microclumping model are not 
satisfied (e.g. if clumps are optically thick for the investigated 
diagnostic), the line formation will depend on more structure 
parameters than just f c \. Thus, relaxing the microclumping ap- 
proximation means that we must consider additional parameters 
when describing the structured wind. These parameters (for a 
two component medium) were defined and discussed in PapeiT, 
and are listed in Table [2] We stress that they are essential for the 
radiative transfer in an inhomogeneous medium, and not merely 
'ad-hoc parameters' used in a fitting procedure. 

In addition to the clumping factor f c \ (or alternatively / v ), 
*ic = Pic/ Pel denotes the density ratio of the interclump (ic) to 
clumped (cl) medium. The time interval St (given in units of the 
wind's dynamic time scale, Z?*/Uoo, and not necessarily constant 
throughout the wind) effectively sets the physical distances be- 
tween clumps, also known as the porosity length h (Owocki et al. 
12004 . which in our geometry is given by h = vpSt. Moreover, 
assuming a smooth underlying field of customary yS-type, vp(r) = 
u<x,(l - bjrf with b set by the assumed velocity at the wind base 
^rnin = v(r =1), this time interval sets the velocity separation be- 
tween the clumps Ad =s vpStdvp/dr (Appendix A). Finally, the ra- 
tio of the clump velocity span Sv (as defined in Fig. 2 in Paper I) 
to this velocity separation (representing a velocity filling factor, 
see Appendix A) largely controls how strongly a perturbed ve- 
locity field affects line formatiorQ. In addition to these basic pa- 
rameters, the velocity u i (or radius r c \) at which clumping is as- 
sumed to start also plays an important role for the line formation. 
Note also that the parameters defining these stochastic winds are 
independent of the physical origin to the inhomogeneities. 

The stochastic models should be distinguished from the 
time-dependent RH simulations. In the latter the structure arises 
naturally from following the time evolution of the wind and 
stems directly from the line-driven instability. Thus, the time- 
averaged structure parameters, as functions of radius, are an out- 
come of these simulations (in contrast to the stochastic mod- 
els, where they are used as fundamental parameters defining the 
structured wind). Nonetheless, the exact wind structure still de- 
pends on the chosen initial conditions, for example on whether 
the instability is self-excited or triggered by some excitation 
mechanism (the latter is done here, see Table [TJ. Finally, as 
shown in Paper I, by choosing a suitable set of structure parame- 
ters one can reconcile spectrum synthesis results stemming from 
the stochastic models with those from RH simulations. 



1 In this paper, we do not consider the 'jump velocity parameter', v s , 
defined in Paper I, since it was shown there that this parameter mainly 
influences the formation of very strong saturated lines, which are not 
considered here. In our applications in Sect. [4] we simply set v-Jvp = 
0.15, which was found to be a prototypical value in Paperl. 



2.2. Radiative transfer 

For resonance lines we use the Monte-Carlo code described in 
Paperl, but a new radiative transfer code has been developed 
for the synthesis of wind recombination lines presented here. 
We investigate the O star recombination lines H a and He n 4686. 
Recall that recombination lines and resonance lines are formed 
differently. First, the optical depths are calculated in different 
ways. For resonance lines, the optical depths may be computed 
via a line-strength parameter, Kq, which is assumed to be con- 
stant throughout the wind and is proportional to the product of 
the mass-loss rate and the abundance (by number), a\ = n jnu, of 
the co nsidered element i. ko may be expressed as (e.g., Pui s et al] 
2008) 



K0 = 



M ne ' /m e c a\ 
R*vl, 4nm H l+4Y He 



(1) 



where e and m e are the electron charge and mass, respectively, c 
the speed of light, ;«h the atomic hydrogen mass, Tee the helium 
abundance, f\ a the transition's oscillator strength, and Ao its rest 
wavelength. The advantage with this definition is that the radial 
Sobolev optical depth in a smooth wind collapses to 



TSob = q 



r 2 udu/dr' 



(2) 



where q is the ionization fraction of the considered ion, and r is 

measured in units of and v in units of Vca- 

For H ff , the analog to kq is the parameter A dPuls et al.ll 19961 
Eqs. 1-3). A is proportional to the mass-loss rate squared and to 
the NLTE departure coefficient, b\, of the lower transition level 
(minus the correction factor for stimulated emission). b\ = tii/n*, 
where n* is the occupation number of level i in LTE with respect 
to the ground state of the next ionization state. In this c ase, the ra- 
dial S obolev optical depth in a smooth wind becomes dPuls et al.l 
1996) 



TSob 



A(r) 



r 4 v 2 dv/dr 



(3) 



In addition to their different optical depths, recombination lines 
are (mainly) formed by recombining ions creating wind photons, 
whereas resonance lines are formed by re-distributing photo- 
spheric stellar continuum radiation by line scattering. Therefore 
the participating atomic levels for recombination lines are rather 
close to LTE with respect to the next ionization state (see Fig.[5]l, 
which means that the departure coefficients and thereby also the 
line source function, S\ oc (e ' kT b\lb u - for these lines 
are basically unaffected by the radiation field and its dilution. In 
turn this al lows us to prescribe the source functions (IPuls et al.l 
1996, l2006h and simply carry out the 'formal integrals' within 
our stochastic and RH winds. In the present work, we assume 
that changes in the NLTE departure coefficients due to optically 
thick clumps can be neglected for recombination-based line for- 
mation, and simply calculate the foj's from NLTE model atmo- 
spheres using the microclumping approximation. Taking the ex- 
ample of Hq, in O stars, this assumption should be reasonable, 
for the H„ departure coefficients in this domain are very close 
to unity and the ionization of hydrogen is complete. However, 
for the case of, e.g., A-supergiants, the assumption no longer 
holds, because in that stellar domain H„'s lower level becomes 
the effective ground state of hydrogen, which means that the 
line transforms to a quasi-res onance line (and thereby that Si 
depends on the radiation field. IPuls et al.lll998h . 
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As described in Paper I, also the ionization fractions q for 
resonance line formation are calculated assuming microclump- 
ing. These fractions are used both within clumps and for the 
interclump medium. The potential feedback effects of optically 
thick clumping on the departure coefficients and ionization frac- 
tions will be investigated by incorporating the analytic methods 
developed in Sect. [3] into suitable NLTE atmosphere codes, and 
reported in a future paper. 

The assumption of prescribed departure coefficients is an 
enormous simplification compared to the UV resonance lines, 
and has enabled us to extend our 2D wind models to 3D when 
modeling recombination line s. In the s ynthes is we follow the 
basic method introduced by iPuls et all d!996l) . which does not 
rely on the Sobolev approximation, with appropriate modifica- 
tions for the line opacities of Hen 4686. A core/halo approach 
is adopted, in which a photospheric profile is used as a lower 
boundary input (at r — 1, with r in units of the stellar radius) 
and the radiative transfer is solved only in the wind. As for the 
resonance lines (Paper I), we assume pure Doppler line broaden- 
ing within the wind, characterized by a thermal speed, u t , given 
throughout the paper in units of i^. 

Our new recombination line code has be en extensively tested 
and showed to yield equivalent results with lPuls et alj d2006l) for 
smooth winds. Also, results based on the microclumping tech- 
nique are reproduced for stochastic as well as RH winds with 
low wind densities, as expected because the clumps then remain 
optically thin. In our applications, we use hydrogen and helium 
occupation num bers calculated by fastwind model atmospheres 
dPuls et al.ll2005l) . under the microclumping approximation, as 
input for the radiative transfer to compute synthetic spectra. 
Photospheric profiles are taken from NLTE calculations of atmo- 
spheres with negligible winds. The consistency between unified 
(meaning a simultaneous treatment of the photosphere and wind) 
model atmosphere calculations and the simplified core/halo ap- 
proach has been verified in the microclumping limit, for recom- 
bination lines as well as for resonance lines. Moreover, we have 
found that averaged recombination line profiles calculated from 
our earlier 2D, stochastic models are almost identical to those 
calculated from our new 3D ones, as was already anticipated for 
the UV resonance lines in Paper I. 

The Hen blend in H„. The star's helium abundance has of 
course been considered in the calculation of the H a wind opac- 
ity, but for simplicity we include the He n blend only in the pho- 
tospheric profile, thus neglecting its direct contribution to the 
wind emission. This results in a slight underestimate of the to- 
tal wind opacity of the line complex. However, by comparing to 
unified model atmosphere calculations that consistently treat the 
He ii blend, we have found that the direct helium contribution is 
low for our typical stars of interest, and in our applications for 
A Cep it can even be neglected. Although sufficient for our pur- 
poses here, this approach should obviously not be generalized, 
because it may yield unrealistic results for stars with parameters 
different from our template star. 

2.3. Geometry 

To c onstruct (pseudo-)3D win ds, we use the 'patch method' 
fromDessart & Owockil d2002l) . A standard right-handed spher- 
ical system (r, ©, <t>) is used, defined relative to a Cartesian set 
(X, Y, Z). However, when computing recombination lines, we no 
longer assume symmetry in the azimuthal (O) direction (as was 
done in Paper I). The lateral scale of coherence in the wind is set 




Fig. 1. Illustration of the wind geometry, see text. A color version 
of this figure is available in the web version. 



by the parameter N@ and by assuming that the physical coher- 
ence lengths in both lateral directions are approximately equal 
(Fig. [TJ. This assumption is reasonable because, within our ap- 
proach, which for example does not include an axis of rotation, 
all observer directions should be alike. Thus, if we desire a co- 
herence scale of 3 degrees, the number of slices in the polar di- 
rection should be N@ - 180/3 = 60 and in the azimuthal di- 
rection N<s, = int[2Af© sin 0], i.e., 2N@ at the equator but fewer 
toward the pole in order to preserve the physical length scales. 
Wind slices are then assigned randomly from a large number of 
spherically symmetric simulations (either RH or stochastic). 

N@ thus enters all our models as an input parameter. Paper I 
showed that this parameter does not change the strengths of 
the resonance lines. More tests have shown that also the ef- 
fects on recombination lines are modest for investigated values. 
Therefore all 3D models in this paper assume N@ = 60, meaning 
a coherence length of 3 degrees at the equator, which is con- 
sistent with observ ational constraints derived from line-profile 
variability analysis dDessart & Ow ocki 2002). Theoretical con- 
straints on N@ are still lacking, and will require a careful treat- 
ment of the late ral radiation transport in RH models. The first 2D 
simulations by iDessart & Owockil d2003l) neglected this trans- 
port and resulted in a laterally fr agmented wind down t o the 
grid scale but the follow-up study (De ssart & Ow ocki 2005b in- 
cluded a simplified 3-ray approach and resulted in larger (but 
un-quantified) lateral coherence scales. 

For recombination line formation, the observer is assumed 
to be located at infinity in the Z u (subscript u denoting a unit 
vector) direction. The geometry is sketched in Fig. Q] We solve 
the radiative transfer using a traditional (P, Z) system for a set of 
P-ray s, each defined by the minimum radial distance to the Z axis 
and by the azimuthal angle <t), which is constant along a given 
ray. If the angle between the ray and the radial coordinate is 6, 
then fi = cos 8 and P — r \J1 - fJ 2 . Thus, for rays in direction Z u 
the radiation angle 6 coincides with the polar coordinate ©, and it 
becomes trivial to calculate the physical locations at which wind- 
slice borders are crossed. The observed flux may then finally 
be computed by performing a double integral of the emergent 
intensity over P and 
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3. Theoretical considerations of resonance and 
recombination line formation in clumpy winds 

Resonance line formation in clumpy hot star winds was dis- 
cussed in detail in Paper I. There we identified an intrinsic cou- 
pling bet ween the effec ts of porosity and vorosity (='velocity 
porosity ', Owocki 2008), which we here further elaborate upon. 
In particular, we propose an analytic formulation of line forma- 
tion in clumped hot star winds (that does not rely on the mi- 
croclumping approximation). As already mentioned in Sect. |2 
the development of such simplified approaches is important for 
properly including effects of optically thick clumping into at- 
mospheric NLTE codes. For recombination lines, we focus on 
H a and discuss impacts from optically thick clumping on its for- 
mation, using our stochastic wind models as well as an extension 
of the analytic treatment developed for the resonance lines. 

3.1. Analytic treatment of resonance lines in clumpy winds 

Throughout this section we assume a smooth velocity field, char- 
acterized by p — 1 (Sect. 12. U . Despite this, the vorosity effect 
will be demonstrated to be important for the line formation (i.e., 
a non-monotonic velocity field is not required for vorosity to 
be at work). Adopting the normalized, dimensionless frequency 
x = (c/VooXAj) - A)/A, and following the basic arguments of 
lOwockil d2008l) . we write the absorption part of a normalized 
resonance line profile, R d x , from a radial ray as (see Appendix 
A) 



/?a,*=6te" TcU + (l-&K 



(4) 



R ax describes the part of the profile that stems from absorption 
of continuum photons released from the photosphere. The to- 
tal line profile is given by R x = R dx + R emyX , where R em , x is 
the re-emission profile. The scattering nature of the resonance 
line source function significantly complicates the formation of 
Rem.x- Therefore, for now we restrict ourselves to discussing R dx 
from a radial ray for these lines. For recombination lines, on the 
other hand, we will lift these restrictions and treat the complete 
line profile (Sect. USV Obviously this must be done also for res- 
onance lines before, e.g., including our analytic formalism in 
NLTE model atmosphere codes (see also Sect. I6.21 >. In any case, 
recall that R dx controls the actual line-profile strengths of reso- 
nance lines, because these are pure scattering lines formed out 
of re-distributed continuum radiation emerging from the photo- 
sphere. 

In Eq.|4]we define £ as the fraction of the velocity field over 
which photons may be absorbed by clumps, with the t's rep- 
resenting the optical depths for the clumped (subscript cl) and 
rarefied (subscript ic) regions. £ describes the essential effects 
of Owocki's vorosity; the first term in Eq. [4] handles the part 
of the line profile emerging from absorption within the clumps, 
whereas the second term handles the part emerging from absorp- 
tion within the interclump medium. What remains then is finding 
an appropriate expression for In Appendix A we argue that a 
reasonable approximation may be 



Ay Av 



fvel(v) + C(v)T](v), 



(5) 



with Av the velocity gap between two clump centers, / ve i the 
velocity filling factor (defined in full analogy with the tradi- 
tional volume filling factor), r\ the effective escape ratio (here 
re-defined from Paper I, see Appendix A), and C a correction 




Fig. 2. Panel a: £ (Eq.|5j as a function of wind velocity. Panel b: 
The relative contributions to £ from / ve i and Cr\. Curves in pan- 
els a and b have been calculated with kq = 5 and thermal speed 
Vt = 0.005 (in units of Panels c andd: Analytic (Eq.[4] solid 
lines) and Monte-Carlo (dashed lines) based absorption part res- 
onance line profiles from clumped winds, as compared to smooth 
results (dotted lines). Clumping starts at v c \ = 0.2 (r c \ = 1.24). 
Panel c: Profiles for three different values of the line strength 
Aft (indicated in the figure), with increased absorption for higher 
values of a:o. Only the ko = 1 profile is not saturated (i.e., does 
not display R d = 0) for smooth models. Panel d: Profiles for 
ato = 5 and different values of u t , as indicated in the figure. No 
Monte-Carlo profile for the v t — case is shown (simply because 
we have not yet developed a 'Sobolev version' of this code). 
Analytic and Monte-Carlo profiles display increased absorption 
for higher values of u t , whereas all smooth profiles have R d = 0. 



factor that depends on the line strength. The condition for inter- 
action for a radial line photon is in the Sobolev approximation 
simply x = v, and, for simplicity, we from now on suppress all 
velocity dependencies of the quantities in Eq. As shown in 
Appendix A, we may write 



/vei — fii 



Sv_ 

6v B 



(6) 



where h is the porosity length of the medium and L = uj(du/dr) 
the (in this case radial) Sobolev length. For the smooth velocity 
field considered in this section, 5u = dvp, which gives / ve i=/v 
Even though the principle effect of the optically thick clumps 
on resonance line formation is a velocity effect governed by / ve i, 
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Eqs. [4]|6] indicate there is also a dependence on spatial porosity 
through the ratio 77 = L/h. This coupling was argued for already 
in Paper I. However, it appears that £ better characterizes the ef- 
fects of clumping in resonance line formation than did our pre- 
vious parametrization (see Appendix A). The optical depths en- 
tering Eq.|4]may be approximated by the corresponding Sobolev 
ones, corrected for the influence of 77 (Appendix A). Then for 
resonance lines, t c i, x ~ T sm , x /£ x and T ic , x = T C K X x ic , where r sm , x 
is given by Eq. [2] Note that all parameters used to define our 
stochastic wind models (Table [2]) enter the expression for R a x , 
illustrating that indeed all these are important for the general 
line formation problem. 

The upper two panels of Fig. [2] plot £ as well as the rela- 
tive contributions from / ve i and Crj for a resonance line with 
line-strength parameter ko = 5. We do not allow for £ values 
above unity (important for the outer wind, see also Appendix 
A), and we also recover the smooth result for velocities lower 
than u c i simply by setting i; % < Vd = 1. All models displayed in 
Fig.|2]were calculated with density structure parameters / c i=4.0, 
(5f=0.5, and jtj c =0.0025. For a smooth model (with ionization 
fraction q = 1, assumed in this section), kq = 5 results in a pro- 
file at the saturation threshold. In the lower two panels we show 
analytic absorption-part line profiles calculated using Eq.|4]and 
profiles calculated using our Monte-Carlo code. To make con- 
sistent comparisons between methods, we accounted only for ra- 
dial photons in the Monte-Carlo simulations. The agreement be- 
tween the methods is very good, lending support to the proposed 
analytic treatment and providing a relatively simple explanation 
for the basic features of the synthetic profiles. 

Evidently, profile-strength reductions can be quite dramatic 
for 'moderately strong' cases such as kq = 5. For the very strong 
kq = 500 line also the interclump medium is optically thick and 
the profiles are therefore saturated (which is a necessity because 
such saturated profiles are observed in hot stars). Note that, if 
6v were much higher than Cv l7 one could neglect the second 
term in Eq. [5] and £ would become independent of the poros- 
ity length. If one also neglects the interclump medium (setting 
x\ c -Q), and assumes that clumps are optically thick throughout 
the entire wind (appropriate for the kq = 5 line), then the ob- 
server in our example would simply receive a constant resid- 
ual flux R d x = 1 - f v = 0.75. Fig. [2}I shows that this gener- 
ally does not hold (even for the idealized case of zero thermal 
speed, the interclump medium still plays a role), demonstrating 
that, along with the velocity filling factor / ve j, in general both 
x; c and r] also help shape the emergent profile for a wide range 
of line strengths and structure parameters. Fig.[2i illustrates the 
importance of accounting for the finite line profile width. Crj 
may not be neglected, even in models with very low, but finite 
thermal velocity, and becomes particularly important toward the 
blue edge of the line profiles. This occurs because the resonance 
zones in the outermost wind become very radially extended. L 
thus grows whereas the distances between the clumps (deter- 
mining h) are unaffected due to the very slowly changing ve- 
locity field. Consequently 77 becomes very high and £ eventually 
reaches unity. Since the smooth kq = 5.0 line is optically thick, 
a 'blue absorption dip' (extensively discussed in Paper I) is cre- 
ated. 

Randomization effects are here neglected because we have 
used a smooth velocity field. When clumps are allowed to have 
velocities higher and lower than those given by the mean ve- 
locity field, overlapping velocity spans of the clumps lead to in- 
creased escape of photospheric photons. The blue absorption dip 
then becomes less prominent than what is displayed in Fig. [2] as 
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Fig. 3. Upper panel: Deviations from the microclumping ap- 
proximation of equivalent widths of synthetic H,, line pro- 
files versus mass-loss rate (see text). Clumping starts at v c \ = 
0.24 and 0.06 (r c i = 1.3 and 1.05), respectively, as indicated 
in the figure. Lower panel: H a line profiles as calculated by 
stochastic, analytic, and microclumped models with fa=9 and 
M=10 x 10~ 6 M Q /yr, and the rest of the stellar and wind param- 
eters as for A Cep. Clumping for all models starts at u c \ = 0.24. 



discussed in Paper I (see also Appendix A, for some comments 
on randomization effects). 

Nevertheless, this section demonstrates that the microclump- 
ing approximation can result in large errors if indeed the wind is 
clumped but the clumps are not optically thin. First applications 
of the analytic formulation are given in Sect. [6] for diagnostics 
of weak wind stars and for the predicted profile-strength ratios 
in resonance line doublets. 

3.2. Recombination lines in clumpy winds 

We now leave the resonance lines behind and turn to the forma- 
tion of recombination lines. We focus on H a , the primary spec- 
troscopic mass-loss diagnostic forO stars. He 11 4686 reacts sim- 
ilarly as H a to clumping in our primary stars of interest (because 
He in is the dominant ion in the line forming regions), and will 
be considered only in our diagnostic study of A Cep (Sect.|4|. 

First we present results from calculating H a line profiles us- 
ing our stochastic 3D wind models. Our main interest is to inves- 
tigate differences with respect to the microclumping model, so 
main results are provided in terms of the deviation of the equiv- 
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alent widths Wa between the two methods, (W^rmc ~ W^/W^rmc, 
as a function of mass-loss rate (here W^mic denotes W,\ as calcu- 
lated from a model assuming microclumping). All models dis- 
cussed in this section were calculated with unity departure co- 
efficients, wind electron and radiation temperatures as for ap- 
proximately A Cep (calibrated using unified NLTE model at- 
mospheres, see lPuls et alJl2006h . and no input photospheric ab- 
sorption profiles. We used structure parameters / c i=9.0, 5t-Q.5, 
Xi c = 0.0025, and a smooth velocity field characterized by (3-1. 

For typical O-supergiants, the equivalent widths of profiles 
calculated from stochastic models are slightly lower than those 
based on the microclumping technique (Fig. [3]). Deviations stem 
from optically thick clumps. The dominating effect is on the 
wind emission of H a photons rather than on the wind absorp- 
tion of photospheric photons (in contrast to resonance lines, see 
previous section). This is because the source function for re- 
combination lines is basically unaffected by the dilution of the 
radiation field, which for relatively strong and hot winds make 
these lines appear in emission and thereby suffer the main effect 
from a clumped wind on the emission part of the line profile. 
Moreover, the p 2 -dependence of recombination line opacity in- 
creases the contrast between the optical depths for the clumps 
and those for the interclump medium, as compared to resonance 
line formation. This lowers the significance of the interclump 
medium and also causes the clump optical depths to decrease 
faster for increasing radii. The latter effect results in clumps that 
are optically thick only in the lower wind regions. Deviations 
from the microclumping limit are therefore more significant for 
cases with earlier onset of clumping. For example, the equiv- 
alent widths for the models with M-2.5 x 10~ 6 M G /yr are re- 
duced by 7 % and 17 % when clumping starts at v c \ = 0.24 and 
v c \ = 0.06, respectively. The effect is thus modest, but notice- 
able. Remember that reductions are measured against models as- 
suming microclumping; the profiles are still much stronger than 
profiles computed from smooth models with the same mass-loss 
rate. 

Our tests show that effects are confined to the line core and 
that the microclumping approximation provides accurate results 
in the line wings. However, Fig. [3] reveals prominent emission 
strength reductions for stronger winds, since then optical depth 
effects become important for ever larger portions of the total 
wind volume. Furthermore, the onset of clumping is irrelevant in 
these strong winds because the majority of the emission emerges 
from radii greater than r c \. This insensitivity to the onset of 
clumping also recovers the scaling invariant for microclumped 
winds (oc -^(/^M, see Sect. 12. 11 1. For typical OB-supergiants, 
however, this scaling does not hold because of the strong opacity 
contrast between wind radii lower than and greater than r c \. Very 
strong recombination lines, such as those displayed in Fig. [3] 
(lower panel), are typically seen in WR stars, for which how- 
ever also a reduced hydrogen content is expected (as well as a 
breakdown of our assumption of an optically thin continuum). 
Nonetheless, our analysis could, of course, be generalized to re- 
combination lines of other chemical species (as has been done 
for He ii 4686 in our application to A Cep), and may point to 
significant optical depth effects in the strong emission peaks 
of stars with very high mass-loss rates. Indeed, lower emis- 
sion peak s in the theoret i cal spectrum of a WR star have been 
found by Oskinova et al.l d2007l) . on the basis of scaling smooth 
opacities using a porosity formalism. However, when deriving 
empirical mass-loss rates from microclumping models of WR 
stars one normally considers also the electro n scattering w ings 
(which are unaffected by microclumping, see iHillierH l 99 1 ). and 
because clumps in these probably are optically thin it may be that 



Table 3. Structure parameters for an empirical stochastic wind 
model of A Cep. 
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lower emission peaks would have a greater effect on the inferred 
clumping factors than on the mass-loss rates. 

Analytic treatment of recombination lines. We can understand 
the reduction in H a emission strengths using the same analytic 
treatment as outlined for resonance lines. Better yet, because the 
source function S is almost unaffected by the radiation field (see 
Sect. 0, we can for recombination lines simulate the total pro- 
file, R x = fl a , x + fl em>x , writing 

/W = - e~ TcU ) + S(l - £0(1 - e~ Ti °*), (7) 

where S is given in units of the continuum intensity and eval- 
uated at the resonance point. R em ^ is much more influenced by 
non-radial photons than is R a x , so accordingly the radial stream- 
ing assumption from the previous section must be relaxed here. 
Details are given in Appendix A. Moreover, the clump and in- 
terclump optical depths now become t c i jX « T srn>x (/ c i/f x ) and 
i"i c ,x = Tci,x*f c > due to the p 2 -dependence of the line opacity, 
with T sm x given by Eq. [3] Already in the previous paragraph we 
mentioned how this lowers the significance of the interclump 
medium in recombination line formation. Actually, tests have 
shown that, in our typical stars of interest, the optical depths in 
the interclump medium are so low that the second term in Eq. [7] 
can safely be neglected. The lower panel of Fig.[3]illustrates that 
profiles computed using the analytic approximation agree very 
well with those computed using our stochastic wind models. 

4. A multi-diagnostic study of A Cep 

We have carried out a detailed study of the Galactic 06 su- 
pergiant A Cep. This star was chosen in part to connect with 
Paper I and in part because it is a well observed and studied ob- 
ject, with significant mass loss, that appears to be less peculiar 
than, e.g., £ Pup. A simultaneous investigation of optical diag- 
nostics and the Pv resonance lines is performed. The ionization 
fractions of Pv and the hydrogen and helium departure coeffi- 
cients (see Fig. [5]l are calculated with the unified model atmo- 
sphere code fastwind, under the microclumping approximation 
and assuming the same (smoothed) clumping factors as in cor- 
responding RH or stochastic models, with stella r and wind pa- 
ramet ers as given in Table Q~| and with a solar ( Asplund et al. 
2005) phosphorus abundance. We account for stellar rotation by 
convolving the emergent synthetic line profiles with a constant 
usin; = 220kms _I , which is the photospheric value. Thus we 
neglect differential rotation in the wind and simply assume that 
stellar rotation only affects the line formation in the wind in such 
a way that we may approximate the same v sin i value as for the 
stellar photosphere (see also Sect. 15. U . Generally, for the line 
profiles studied here, the influence of rotation is important for 
the recombination lines, but not for the resonance l ines. We use 
observed uv fuse spectra fromlFullerton et al.l d2006l) . and optical 
spectra from lMarkova et all (120051) and A. Herrero (described in 
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Fig. 4. Radial clump optical depths for PV, and radial and tan- 
gential ones for H„, as functions of wind velocity for our empir- 
ical, stochastic wind model of A Cep (Table[3]l. 
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Fig. 5. Clumping factors (upper left), Pv ionization fractions (up- 
per right), and He n4686 and H a departure coefficients bi (lower 
left and lower right, respectively) used in the RH and stochastic 
models of A Cep. Mean wind velocities are on all abscissas. 



lHerreroetai1l2000t) . In addition to H„, Hen 4686, and PV, we 
also consider the wind sensitive cores of H^ and H r . However, 
for these diagnostics we rely entirely on the microclumping ap- 
proximation, which because of their low wind optical depths 
should be sufficient. 

4. 1 . Clump optical depths 

The clump optical depth in the wind is the primary quantity 
governing the validity of the microclumping approximation. In 
Appendix A (see also the previous section) we have provided 
estimates of the clumps' radial Sobolev optical depths in reso- 
nance and recombination lines. However, in our stochastic mod- 
els, clumps do not always cover a complete resonance zone, so 
the Sobolev optical depths must be replaced by optical depth cal- 
culations including the actual line profile (cf. Paper I, Eq. A. 12, 
and note also that these 'exact' optical depth calculations then 
involve / ve i rather than see Appendix A). Within our stochas- 
tic wind models, the radial extent of a clump is l r = vpdtf v , and 
therefore, by transforming to the corresponding velocity width, 
we may readily calculate the 'actual' clump optical depth t q \. 

Fig.|4]shows t c i for Pv and H a , as calculated from our empir- 
ical, stochastic model of A Cep (TableO. The figure shows that 
the radial r c i is significantly higher for Pv than for H ff and, more- 
over, that the linear dependence on the density for resonance 
lines (as opposed to the quadratic dependence of recombination 
lines) causes clumps to remain optically thick in Pv throughout 
almost the entire wind. The figure also illustrates how tangential 
H„ photons have higher optical depths than radial ones in the 
line forming regions. Since non-radial photons are important for 
recombination lines (Sect. 13.2b . this enhances the effects from 
optically thick clumping on the H„ line formation in our empir- 
ical A Cep model (as seen in Fig. [7ji. The differences in clump 
optical depths for radial (p — 1) and tangential (p = 0) photons 
stem from the dependence t c i oc [dv/drp 2 +v/r(l-p 2 )]~ l , and the 
fact that v/r is significantly lower than du/dr in the relevant wind 
regions. In any case, based on these simple estimates, one might 
expect that the basic results of Sect. [3] should hold in diagnostic 
applications of typical O stars. That is, H a should be affected by 



optically thick clumping only in the line core, whereas resonance 
lines should be much more affected over the entire line profile. 

4.2. Constraints from inhomogeneous 
radiation-hydrodynamic models 

Fig.|6]displays line profiles calculated from our RH model of A 
Cep. Consistent fits of the observed diagnostics are not achieved. 
The H„ line wings are reasonably well reproduced but the core 
emission is much too low. The Pv profiles are, actually, better re- 
produced, although stronger than observed toward the blue edge 
of the line complex (the 'blue edge absorption dip' problem, 
see Sect. 13. U . The reasonable Pv fits are due both to adopt- 
ing a rather low mass-loss rate for A Cep (see Table [TJ and to 
lower velocity spans in these RH models than in those analyzed 
in Paper The mass-loss rate was essentially chosen from a 
best compromise when considering the complete diagnostic set, 
however note that a consistent fit to all diagnostics could not be 
achieved, independent of which mass-loss rate was adopted, as 
now discussed. 

The apparent mismatch between H a emission in the core and 
in the wings occurs because f c \ increases rather slowly with in- 
creasing velocity (Fig. 0, which for a given mass-loss rate im- 
plies that the optical depths in the H a core forming regions are 
too low as compared to the optical depths in the wing forming re- 
gions. He ii 4686 is subject to the same mismatch as H a , and also 
the cores of H^ and H y are deeper than observed. The latter fea- 
ture occurs because the photospheric absorption profiles are not 
sufficiently re-filled by emission from the only weakly clumped 
inner wind. Thus, the optical wind diagnostics all indicate that 
the clumping factor as a function of velocity in A Cep differs 
from that predicted by t he RH simulations (see also iPuls et af] 
2006; lBouret et al.1 l2008). On the other hand, any significant in- 
crease in the mass-loss rate to obtain a better fit of the higher 
Balmer lines and the core of H a would produce stronger than 
observed H a and Hen 4686 line wings (as illustrated for H a in 
Fig.[6]l and, vice versa, a reduction of the mass-loss rate to obtain 

2 The exact reasons for the lower spans are still under investigation, 
and will be reported in a future paper. 
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Fig. 6. Observed and synthetic line profiles for A Cep. Dotted lines are the observations. Solid line profiles are calculated from 
the inhomogeneous radiation-hydrodynamic model of A Cep (Table [TJ, and dashed lines from a corresponding fastwind model 
including microclumping. The long-dashed line in the upper left panel is from a RH model in which the density has been scaled to 
mimic an increase in the mass-loss rate by 50 %. 



a better fit of the (blue edge of the) PV lines would produce too 
weak wings. 

Comparison with the microclumping technique. We now com- 
pare results from above with those from a microclumped fast- 
wind model having the same (smoothed) clumping factors as the 
RH model. The Pv profiles calculated using the fastwind model 
are stronger than those calculated using the RH model. We may 
characterize this difference by the difference in the equivalent 
widths Wa of the absorption parts of t he profiles. Wa is roughly 
15 % lower for the RH model (see also IOwockil20 08) . However, 
this moderate reduction in profile strength actually corresponds 
to a reduction in the mass-loss rate by a factor of approximately 
two, because of the resonance lines' slow response to mass loss. 

Resonance line profiles stemming from the RH and micro- 
clumping models also display different line shapes. For RH 
models, significant velocity overlaps stemming from the non- 
monotonic velocity field ensure that the observed flux at the 
blue side of the line center is accurately reproduced without in- 
voking any artificial and highly supersonic 'microturbulence', 



as must be done when using smooth as well as microclumping 
wind models. Although not analyzed here, also the absorption 
at velocities >v x of saturated resonance lines may be repro- 
duced by_J3Jimc<lels without invoking additional microturbu- 
lence dPuls et al.ll 19931 Paper I). For H a and Hen 4686, the RH 
and microclumping models yield almost identical results. This 
occurs because clumps are optically thin in these diagnostics 
throughout almost the entire wind, due to the slow increase of 
/ci with mean wind velocity, which in turn results in wind den- 
sities in the inner wind unable to produce optically thick clumps 
(compare to the empirical models in the following section). 



4.3. Constraints from empirical stochastic models 

Clearly, the RH models fail to deliver satisfactory line profiles 
when their structures are confronted with UV and optical wind 
diagnostics in parallel. Here we use our stochastic models to 
modify the wind-structure parameters and show how the results 
then may be reconciled. This is a first attempt toward our long- 
term aim of using consistent multi-diagnostic studies to obtain 
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Fig. 7. As Fig.|6l but using our stochastic models (solids) with corresponding inferred empirical structure parameters (see Table [3] 
and text). The assumed mass-loss rate is the same as for the RH model of A Cep, see TableQ] 



unique views of empirical mass-loss rates and structure proper- 
ties of hot star winds. 

The mass-loss rate is determined by a best fit to the com- 
plete diagnostic set (giving highest weight to the optical hydro- 
gen lines). We derive the same rate as was previously adopted 
for the RH model of A Cep. This rate (M=1.5 x 1(T 6 M /yr) is 
approximately two times lower than the corre sponding fheoreti- 
cal one obtained using the mass-loss recipe in lVink et al.l (12000) 
(M=3.2 x lO~ 6 M /yr). In the outermost wind, we for now ad- 
here to the con straints on f c \ derived from radio emission by 
iPuls et all (120061) . scaled with respect to the mass-loss rate de- 
rived here. In the inner wind, both the distinct shape of H ( , in 
A Cep3 and the cores of the higher Balmer lines may be used as 
tracers of structure. The H ff absorption trough followed by the 
steep incline to rather strong emission can only be reproduced 
by our models if clumping is assumed to start quite late (see 
also lPuls et al.l 2006: Bou ret et al.l2 008). at a velocity marginally 
lower than predicted by the RH models, however with a much 
steeper increase with velocity (see Fig.|5]l. Also, in the particular 
case of A Cep, the upper limit of the mass-loss rate derived by 



3 which only resembles the P Cygni shapes of the UV resonance 
lines, since it is formed differently. 



IPuls et alJ d2006l) (M=3.0 X 10~ 6 M /yr, inferred by assuming a 
smooth outermost radio emitting wind) results in densities so 
high in the lowermost wind that the H a trough never reaches be- 
low the continuum flux. Moreover, additional constraints come 
from the cores of the higher Balmer lines; the higher the den- 
sities in the lowermost wind, the stronger the re-filling of the 
photospheric absor ption profil e by wind emission. Here as well 
the upper limit from lPuls et al.l provides shallower than observed 
line cores. Thus, if we require f c i=l at the wind base, and if our 
interpretation of the abrupt shift from absorption to emission in 
H a as due to clumping is correct, rather tight constraints on the 
mass-loss rate may be obtained using only optical diagnostics. 

The Hq, time series of Marko va et al.l ((2005) reveal that both 
the height of the emission peak and the depth of the absorp- 
tion trough depend on the observational snapshot, variations can 
reach 0.04 in residual flux units. Therefore it is not critical that 
neither the peak nor the trough is perfectly reproduced by our 
models in Fig|7] (which displays a 'representative' observational 
snapshot). On the other hand, the observations do not indicate 
any significant variation in the position of the emission peak. 
This might be an issue, because the late onset of clumping red- 
shifts the emission peak too much (at least when neglecting 
differential rotation, see Sect. 15. U . whereas an earlier onset of 
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clumping fails to produce an absorption trough. The offset in the 
position of the emission peak is larger than the estimated uncer- 
tainty in the radial velocity correction, which may indicate that 
clumping is only partly responsible for the shape of the H a core. 
Indeed, other interpretations have been suggested, and we com- 
ment on this in Sect. 15.11 

The line shape of Hen 4686 is well reproduced by our 
stochastic models, but not the emission strength. The line re- 
acts similarly to clumping as H a . In order to increase the central 
emission to the observed level we would have to raise the clump- 
ing factor in the inner wind even more, which in turn would pro- 
duce stronger than observed H a emission as well as shallower 
than observed H^ and H y cores. Since hydrogen generally has 
more reliable and robust departure coefficients than helium, we 
have given higher weights to fits of hydrogen lines. Interestingly, 
He II 4686 shows a similar offset as H„ in the position of the 
emission peak. 

The PV resonance lines are much more sensitive to the wind 
structure parameters (see Sect. 13.11 ) than to the mass-loss rate. 
Hence these lines should be used only as a consistency check of 
mass-loss rates derived from other diagnostics. Using the struc- 
ture parameters given in Table |2l our stochastic models yield 
reasonable fits of the PV lines. We use values of 5t and x; c as 
in Paper I, including a higher 6t in the outer wind to account for 
clump-clump collisions, but are able to adopt a higher value of 
\Sv/6vp\, which however is still lower than predicted by the RH 
models. This higher value stems from that we here consider also 
optical diagnostics and from these derive a lower mass-loss rate 
and higher clumping factors than what was assumed in Paper I, 
essentially meaning that larger velocity spans then can be used 
when fitting the Pv lines. 

fa in the inner wind is drastically different from that pre- 
dicted by our RH model for A Cep (FigO, and indicates that 
present-day RH simulations fail to predict observationally in- 
ferred clumping factors, at least for the inner wind. Regarding 
the outermost wind, let us point out that the RH simulations 
used here only exte nd to r » 35, at which f c i is still decreas- 
ing. Simulations by iRunacres & Owockil d2002l) . which extend 
to much larger radii, indicate that the clumping factor settles at 
« 4 in the outermost wind. f c \ a; 4 is consistent with our derived 
m ass-loss rate and t he constraints from radio emission derived 
bv iPuls et alJ (120061) (see above). This suggests that the outer- 
most wind is better simulated by current RH models than the 
inner one. 



Comparison with the microclumping technique. Here we 
compare the stochastic models from above with microclumped 
models calculated with the same clumping factors. When us- 
ing the microclumping technique, the PV resonance lines are 
not directly affected by the structured wind. The mass-loss rate 
adopted in the previous paragraph then produces much too 
strong absorption in these lines, see Fig. [7] Moreover, the high 
clumping factor in the inner wind adopted in our stochastic mod- 
els results in so high densities that the clumps become optically 
thick in H„ and He n 4686 as well. This generally leads to weaker 
emission for the stochastic models than for the microclumped 
ones (Sect. 13.21 ). and / c i's drastic increase from 1 to 28 makes 
the deviation from the microclumping approximation prominent 
in this particular case. We have confirmed that the same emission 
strength reduction results when using our simplified analytic ap- 
proach (Sect. [3~!2b . which supports the rather strong emission re- 
duction that we find in the H,, core as well as indicates that our 



analytic approach indeed might be a promising tool for a consis- 
tent implementation into atmospheric NLTE codes. 

In order to obtain reasonable fits of the PV lines within the 
microclumping approximation we had to lower the mass-loss 
rate significantly, to M-0.4x l(T 6 M /yr (this is the so-called 
'Pv problem', see also Full erton et alj 2006). In turn this meant 
that extreme clumping factors, f c \ ~ 400, in the inner wind were 
required to meet the observed amount of H a wind emission. 
However, we have not been able to achieve a consistent fit of the 
optical diagnostics using these highly microclumped fastwind 
models: if for example H a is fitted then the He n4686 emission 
is much too weak. Overall, the results in this section support 
the view that the extremely low empirical mass-loss rates previ- 
ously indicated from Pv might be a consequence of neglecting 
optically think clumping when synthesizing resonance lines. 

5. Discussion 

5. 1. Are O star mass-loss rates reliable? 

Theoretical rates. The time/spatial averaged mass-loss rate of 
our A Cep RH model differs from the rate of the correspond- 
ing smooth start model (used for initialization) by less than 5 %. 
From this one might expect that the clumped stellar wind should 
not significantly affect theoretical mass-loss rat es based on the 
line-driven wind theory. However, RTti cka et al.l d2008l) (see also 
Muijres et al. 2010, submitted to A&A) made some first tests 
and included wind inhomogeneities in a (steady-state) theoreti- 
cal wind model of an O star. They found that the predicted mass- 
loss rate increased when clumps were assumed to be optically 
thin, because of increased recombination rates that shifted the 
ionization balance to lower ionic states with more effective driv- 
ing lines. On the other hand, their tentative attempts to account 
for optically thick clumps in the continuum opacity as well as 
for clumps with longer length scales than the Sobolev length re- 
duced the line force and led to lower predicted rates. 

The reduced profile strengths of resonance lines (which are 
the main drivers of the wind) found here should in principle also 
reduce the line driving in theoretical steady-state wind models, 
but let us point out that many lines that significantly contribute to 
the total driving force might still be saturated because of the non- 
void interclump medium. Nevertheless, it is clear that a thorough 
investigation of the impact of clumping on predicted mass-loss 
rates is urgently needed. The mass-loss rate for A Cep derived 
here is approximately a factor of two lo wer than the theor etical 
rate predicted by the mass-loss recipe in lVink et alJ yOOO). 

Empirical rates. Our empirical mass-loss rate for A Cep is 4.5 
times lower than the rate inferred from sy nthesizing H„ using a 
smooth wind model (Repolus t et al.ll2004l) . The best constraints 
on the mass-loss rate in our analysis come from the distinct 
shape of the H„ line core and the higher Balmer lines (Sect. l4~3l l. 
Rotation in our models is treated by the standard convolution 
procedure. But A Cep is a fast rotator (Table [TJ, so differential 
rotation might influence the formation of the lin e profiles, par- 
ticula rly the Hq, core. Bouret et al. (published in Bresolin et al. 
2008) found that the H„ line in £ Pup can be fitted by assum- 
ing that clumping starts close to the wind base, ;/ differential 
rotation is treated consistently. Since f Pup and A Cep display 
similar H a profiles, it is possible that the same effect could be at 
work also in the latter star, and thereby that the rather late onset 
of and the rapid increase of clumping in our stochastic model of 
A Cep could be somewhat exaggerated. Naturally, this could then 
also affect the inferred mass-loss rate, since with a modified run 
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of the clumping factor another rate might be required to obtain a 
simultaneous fit of the observed diagnostic lines. 

The influence of X- ray and xuv/euv radiation as created by 
shocked wind regions (Feldmeie r et al.l fl997l) on the occupa- 
tion numbers is not included in our analysis. These contributors 
are n ot important for calcu lations of hydrogen occupation num- 
bers (Pauldrach et a l.l2001l) . but their significa nce for the ioniza- 
tion fractions of phosphorus is st ill debated (Rrticka & Kubat 
2009; IWaldron & Cassineifl l2010h . We have used the alterna- 
tive unified atmospheric code WM-Basic (lPauldrach et al. 
which treats X-ray and xuv/euv radiation but not wind clump- 
ing, to estimate the impact of X-rays on the Pv ionization frac- 
tions. We find that effects are negligible at wind velocities lower 
than v/vca « 0.5 but profound at higher velocities, with the 
Pv ionization fraction significantly reduced when X-rays (and 
of course the corresponding xuv/euv radiation tail) are included. 
This suggests that a proper treatment of these hot radiation bands 
might resolve the earlier discussed 'blue absorption dip' prob- 
lem, which is clearly visible in the Pv line profiles calculated 
from RH models (Fig. [6] but note that we overcame this prob- 
lem in our stochastic models by increasing the distances between 
clumps in the outermost wind, see TableO. 

5.2. Structure properties of the clumped wind 

We identify two main problems when confronting synthetic 
spectra from the time-dependent RH simulations of the line- 
driven instability with observed lines in the UV and optical: i) 
the absorption toward the blue edge of unsaturated UV reso- 
nance lines is too deep in the simulations, and ii) the emission 
in the core of H a is much too weak as compared to the emission 
in the wings. The first problem is related to the high predicted 
velocity spans in the RH models, and was extensively discussed 
already in Paper! Moreover, in Sect. l5.ll we commented on that 
even if the large velocity spans turn out to be stable features, this 
problem might be overcome by a proper treatment of X-rays in 
the calculations of ionization fractions. 

The second problem arises because the predicted clumping 
factors in the inner wind are too low as compared to those in 
the outer wind (Fig. |5j. However, let us point out that veloc- 
ity as well as density perturbations in the inner wind of our 
RH simulation may be overly damped, because we use the so- 
called smooth source function (SSF) approximation when cal- 
culating the contribution to the line force from the diffuse, scat- 
tered radiation field. In simulations that relax the SSF approxi- 
mation and account for gradients in the perturbed source func- 
tion (via an 'escape-integral source function' formulation, EISF, 
lOwocki & Pulslfl996l 1 1999b . the structure in the inner wind is 
more pronounced and also develops closer to the photosphere. 

In any case, however, it is questionable if self-excited in- 
stability simulations will be able to reproduce the observed 
clumping patterns (which have been found also in earlier in- 
vestigations based on the microc lumping approximation, e.g., 
Bouret et al. 2 0051; iPuls et af]|2006l) . especially considering that 
our RH model of A Cep actually already is triggered (Table []]), 
using Langevin perturbation s mimicking photospheric turbu- 
lence (Feldm eier et alj [l997). Thus, while observations tracing 
the outer wind seem to confirm the structures predicted by the 
line-driven instability, observations tracing the inner wind might 
require the consideration of an additional triggering mechanism 
to be reproduced, which perhaps must be stronge r than what 
is currently assumed. For example, ICantiello et al.l (|2009) pro- 
posed that gravity and/or acoustic waves emitted in sub-surface 
convection zones may travel through the radiative layer and 




<q>dM/dt [M Sun /yr] 

Fig. 8. Equivalent widths, Wa (normalized to the value for a sat- 
urated line), for the absorption part of the Nv resonance line at 
1240 A, as functions of the product of the ionization fraction of 
Nv, (q), and the mass-loss rate. The solid line is calculated from 
smooth models and the dashed line from structured ones. The 
black dots denote W/s for models corresponding to a smooth 
model with (q)M - 10~ 9 , see text. 

induce clumping already at the wind base. However, regard- 
ing gravity waves, it is not certain that these would have high 
enough frequencies (i.e., higher than the atmosphere's acoustic 
cutoff frequency) that they can be radially transported through 
the wind. Another possibility for a strong clumping trigger might 
be non-radial pulsations in the photosphere. Certainly it would 
be valuable to investigate to what extent such triggers, within 
a line-driven instability simulation using the EISF formulation, 
could produce clumping patterns in the inner wind more com- 
patible with the observations. 

6. Additional considerations 

In this section, we discuss two applications for the analytic 
formulation of line formation in clumpy winds presented in 
Sect. IP 

6. 1 . Weak wind stars 

The so-called weak wind problem is associated with observa- 
tions of (primarily) O-dwarfs of late types, which appear to have 
mass-loss rates much lower than what is predicted by the line- 
driven wind theory, and also much lower than other 'normal' O 
stars of earlier spectral types. However, a major problem with 
wind diagnostics in this domain is that the primary optical di- 
agnostic, H a , becomes insensitive to changes in the mass-loss 
rates, so that only upper limits can be inferred from this line. 
Therefore one must for these objects quite often rely solely on 
the intrinsically stronger UV resonance lines. For a comprehen - 
sive discussion on the weak wind problem, see Pu is et al. (2008). 

In the following, we demonstrate the potential impact of op- 
tically thick clumping on diagnostic resonance lines in weak 
wind stars using the analytic formulation developed in Sect. I3.ll 
We use one component of t he Nv doublet at 124 A, assume 
a solar nitrogen abundance dAsplund et alj 120051) . and take a 
generic O-dwarf with parameters /?*=8.0/? Q , Uo^lSOOkms -1 , 
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and Yn e =0.1. Since in this section we only discuss predictions 
for the product of mass-loss rate and ionization fraction q for 
resonance lines, no effective temperature needs to be specified 
(see Eqs. [T E1>. The Nv doublet wa s among the lines utilized in 
the study of Marcol ino et al.l (120091) . and also our chosen param- 
eters correspond well to the parameters for the five stars ana- 
lyzed and found to have very weak winds (more than an order 
of magnitude lower than predicted by theory) in that study. To 
avoid problems regarding the onset of clumping and the afore- 
mentioned 'blue absorption dip', we consider only the veloc- 
ity interval v/voo = 0.25 - 0.75. Absorption-part line profiles 
for structured winds are calculated using Eq.|4]and adopting the 
same structure parameters as in Sect. 13. IK £,=0.25. Xi c =0.0025, 
(5f=0.5, and a smooth '/3=V velocity field). 

Fig.|8]shows the curve-of-growths for structured and smooth 
models, respectively, as functions of the mean ionization frac- 
tion of Nv times the mass-loss rate, (q)M. Clearly, mass-loss 
rates derived from smooth models may be severely underesti- 
mated also for stars with weak winds. For example, if we for 
this star were to infer {q)M = 10~ 9 M /yr from a smooth model, 
the corresponding rate inferred from a structured one would be 
«<7>M) struc = 3.8 x 10- 8 M /yr = 38((q)M) smooth (see Fig.|8j. 
Thus, if using smooth models (or microclumped, since micro- 
clumping has no effect on the resonance lines as long as no sig- 
nificant changes occur in the ionization fractions), one could eas- 
ily derive mass-loss rates more than an order of magnitude lower 
than corresponding rates derived from structured models (see 
also [Oskinova et al. 2007 and Paper I), and thereby one could 
also misinterpret observations as suggesting that mass-loss rates 
are much lower than predicted by theory. 

We emphasize, however, that this simple example merely 
demonstrates how optically thick clumping might be impor- 
tant also for resonance line diagnostics in so-called weak wind 
stars, and that, the winds are clumped, one must be care- 
ful not to simply assume that strongly de-saturated resonance 
lines also imply optically thin clumps. The actual mass-loss 
reductions will depend critically both on the assumed ioniza- 
tion fractions and on the adopted structure parameters. Thus, a 
multi-diagnostic study (to constrain the structure parameters), 
including a detailed consideration of X-rays (to obtain reliable 
ionization fractions), is required for more quantitative results. 
Nevertheless, we may safely say that, because of these inher- 
ent problems in UV line diagnostics, it is important to put fur- 
ther constraints on the weak wind problem by exploiting other 
diagnostics that are sensitive to mass loss but neither have opti- 
cally thick clumps nor are affected by X-rays (as is probably true 
for, e.g., the in frared Br g line, Najarro et al. 2010, submitted to 
A&A, see also lPuls et aD2009l) . 



6.2. Resonance line doublets 



iMassa et~ai1 (|2008) pointed out that additional empirical con- 
straints on wind structure may be obtained by considering the 
observed profile-strength ratios of resonance line doublets. The 
line strength parameter, kq, of such doublets is in proportion to 
the oscillator strengths of the individual components, /, which 
for the cases of interest here are / b // r = 2, with superscripts 
b and r denoting the blue and red line components, respec- 
tively. However, if clumps are optically thick for the investigated 
lines, the resulting profile-strength ratio may deviate quite sig- 
nificantly from the one implied by smooth modeling (see dis- 
cussion in Paper I). For example, in the case of very optically 
thick clumps and a void interclump medium, Eq. [4] simply gives 
R a> x = 1 - £ x , i.e. the profile strength becomes independent of 



kq. The analogy for continuum diagnostics, or for line diagnos- 
tics in a non-accelerating medium, is the well-known result that 
for a medium consisting of infinitely dense absorbers embedded 
in a vacuum, the effective opacity is independent of the atomic 
opacity (see footnote 4 in Appendix A). Also for such a situation 
would the inferred profile-strength ratio be exactly one. 

A major advantage of this line diagnostic is that the depen- 
dence on X-rays sh ould cancel out. Recently, iPrinia & M assa 
(2010) extended the IMassa et"aT1 work to include a large num- 
ber of B supergiants, for which they, from the Si iv /i/U400 res- 
onance doublet, derived empirical line-strength ratios, us- 
ing smooth wind models. The stars showed a wide spread be- 
tween unity and the predicted factor of two, with the major- 
ity of them lying in the range 1.0 to 1.5, and with an overall 
mean of 1.46 (standard deviation ~0.31). In the following, we 
shall discuss this diagnostic under the assumption that the dou- 
blet components are well separated, so that each component can 
be treated as a single line, which is reasonable for, e.g., the just 
mentioned silicon lines in typical B-supergiants and for Pv in 
OB-stars. 

We now show that our analytic formulation for resonance 
line formation indeed predicts profile-strength ratios on the same 
order as those discussed above. Following the preceding section, 
we assume a solar abundance for silicon, make use of a generic 
B-supergiant with fl*=3O.O/? , ^=800 km s" 1 , and Y He = 0.1, 
adopt the same structure parameters as in the previous section, 
and consider only the velocity interval v/v^ = 0.25 - 0.75. We 
then assume that for this generic star we derive (q)M = 5 x 10~ 9 
from the Siiv resonance doublet formed in a structured wind 
model. By once more exploiting the curve-of-growth (as in 
Fig. [8] but now for the two components of Si iv), we can then 
easily translate the structured results to corresponding smooth 
ones. We find a ratio (ftn/*^)smooth ~ 1 - 4, whi ch agrees well with 
the results derived by Prinja & Massa| d2010h . 

The doublet ratios are, in fact, almost ideal diagnostics re- 
garding structure properties, since all other dependencies cancel 
out. Therefore ratios deviating from two might be the cleanest 
indirect signatures of optically thick clumping that we presently 
have, and may in principle be used to extract empirical informa- 
tion on the behavior of We write the ratio of the blue and red 
absorption-part line profile at frequency x as 



(8) 



Generally, this equation can be solved for £ x only if the line op- 
tical depths and the interclump densities are known (the latter 
for example from observations of saturated resonance lines, see 
Paper I). However, under certain circumstances we can eliminate 
the need for such external knowledge. For example, assuming 
that all clumps are optically thick, we may write 



& = 1- 



(* f a,x) 2 

R b ' 



(9) 



Applying the last expression to our line profiles computed for 
Si iv using Eq.|4]reveals a mean value of £ = 0.48 in a velocity 
bin v/vco = 0.4 - 0.5, which agrees well with the actual mean 
(calculated from the assumed structure parameters), ^ = 0.51. 
Thus, this approximation can provide a quite good direct empir- 
ical mapping of ^, without any knowledge about optical depths 
etc. Another case for which the profile-strength ratio can be di- 
rectly related to £ is that of a completely transparent background 
medium (i.e., a void interclump medium in our case). That lim- 
iting case of Eq. [8] has been long recognized and used by the 
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quasar community (e.g jGangulv et aOl999l) . for the formation 
of intrinsic, narrow absorption-line doublets. 

However, let us point out that this theoretical example only 
demonstrates that our basic formalism appears reasonable. In a 
real application, there will be a contribution also from the re- 
emission part of the line profile, i.e., what we actually measure 
from an observation is the total line profile R x = >? ax + 7? em x . 
Thus, to empirically infer £ x from Eq.|9] (which involves R a x = 
Rx - fi e m,x), we must either simply neglect the re-emission con- 
tribution (which generally will not be possible) or actually cal- 
culate R em ,x, as predicted by a structured wind model. For reso- 
nance lines (as opposed to recombination lines, see Sect. 13.21 1. a 
simplified approach for /? e m,x in clumpy winds is still to be de- 
veloped; it is a very demanding task because of the source func- 
tion's scattering nature. In principle though, a treatment corre- 
sponding to the 'smooth source function' formalism used in our 
time-dependent RH simulations (see Sect. 15.2b might be a rea- 
sonable first approximation. 

7. Summary and future work 

We investigate diagnostic features for deriving mass-loss rates 
from the clumped winds of hot, massive stars, without relying 
on the microclumping approximation. It is found that present- 
day RH simulations of the line-driven instability are not able to 
consistently fit the UV and optical diagnostics in a prototypi- 
cal O-supergiant. By creating empirical stochastic wind models, 
we achieve consistent fits mainly by increasing the clumping 
in the inner wind. A mass-loss rate is derived that is approxi- 
mately a factor of two lower than what is predicted by theory. 
The best constraints come from the optical diagnostics. The UV 
resonance lines are much more sensitive to the wind's structure 
parameters (i.e., to the clumping factor, the interclump medium 
density, etc.) than to the mass-loss rate, and should, thus, not be 
the preferred choice when deriving empirical mass-loss rates. 

We discuss both recombination line and resonance line for- 
mation in detail. Resonance lines always suffer the effects of 
optically thick clumping in typical diagnostic lines, and their 
profiles are thereby weaker for models with a detailed treat- 
ment of clumping than for models that rely on the microclump- 
ing approximation. Recombination lines are less affected be- 
cause of the lower optical depths in typical diagnostic lines. 
However, emission strength reductions as compared to micro- 
clumped models are significant for stars with high mass-loss 
rates (e.g., Wolf-Rayet stars) and can be so for O stars as well, 
for example, strong clumping is present in the lower wind, as 
illustrated by our diagnostic study of A Cep. 

An analytic method to model these lines in clumpy winds, 
without any restriction to microclumping, is suggested and 
shown to yield results consistent with those from detailed 
stochastic models. Some first results are given, illustrating the 
potential significance of optically thick clumps for diagnostic 
lines in weak wind stars, and confirming recent results that 
profile-strength ratios of resonance line doublets may be used 
as tracers of wind structure and optically thick clumping. We in- 
tend to refine this method and incorporate it into suitable NLTE 
unified atmospheric codes, in order to investigate effects of opti- 
cally thick clumping on the occupation numbers. 

It is pivotal that 3D, time-dependent RH models of the line- 
driven instability be developed, with an adequate treatment of 
the 3D radiation transport. New models are required to investi- 
gate whether the structure predicted by present-day simulations 
is stable or a consequence of current physical assumptions and 
simplifications. 
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Appendix A: Analytic treatment of line formation in 
clumped hot star winds 

Resonance lines. We propose to write the absorption part of a 
resonance line formed (from a radial ray) in a clumped wind as 



(A.l) 



where £ is defined as the fraction of the velocity field over which 
photons may be absorbed by clumps, the optical depths are those 
for the clumped (subscript cl) and rarefied (subscript ic) medium, 
and dependencies on the normalized, dimensionless frequency x 
have been sup pressed for simp licity (cf. Sect. l3.ll in main paper). 

Following Owockil d2008l) we define the velocity filling fac- 
tor / ve i as the fraction of the velocity field covered by clumps 
(in full analogy with the volume filling factor / v ). That is, f, e \ is 
the ratio of the velocity span of the clump, 5v, to the velocity 
separation between two clump centers, Au, 



/vel = Ay" 



(A.2) 



In our stochastic models we have the clump velocity span 6v » 
\5v 1 5vp\{dvp I dr)6r and from the definition of / v (see Paper I) 6r a 
/ v Ar, with Ar = vp6t the radial distance between two clump 
centers. Similarly one may approximate Ay « (dvp/dr)Ar, which 
leads to 



8v 

Jvel ~ I — l/v 
OVr 



(A3) 



simply modify £ in order to account for the essential effects. We 
write 



£ 



Sv + Cv t 
Av 



(A.4) 



where the factor Cv l now represents a sort of correction to the 
limiting case of 6v » Cv v A linear addition is chosen because 
the basic equation determining whether or not a photon actually 
can be absorbed (i.e., whether or not it is located within its reso- 
nance zone) is x cm f = x b s - v, with jc cm f and x b. s the co-moving 
and observer's frame frequencies, respectively. 

The factor C accounts for the fact that the 'effective reso- 
nance zone' over which photons can be absorbed by clumps is 
larger than that provided by u t (at least for relatively strong lines). 
Photon absorption at x within clumps is given by the distribution 
function e~ Tclx , with expectation value t c i jX = 1. Therefore we 
may estimate C using the 'effective profile width', determined 
by solving for the co-moving frame frequency at which unity 
optical depth is reached, // a clump is present, 



1 - erf[x cmf /u t ] 
Td ™ = 1, 



(A.5) 



where erf is the error function. The effective profile width then 
is C — 2x cm flv t , where x cm f is given by the solution to Eq. IA.5l 
Note that C now is allowed to be velocity dependent, C — > C(v). 
In addition, the expression for the clump optical depth should 
now be modified, T sm // ve i — » r sm /£ to account for the fact that 
individual clumps no longer cover a complete resonance zone. 
With C determined we can cast £ in the convenient form 



Thus, a smooth velocity law (6v—6vp) implies f ve i=f v . The ab- 
solute value sign becomes important in any compressive wind 
region, since optical depths must always be positive (see below). 

Actually, Eq. I A. II is in form equivalent to the analytic trans- 
fer solution derived bv lLevermore et aD dl986l) . for the ensemble 
averaged intensity in a two-phase [i - A, B] Markovian model of 
a static purely continuum absorbing medium in the limit that the 
characteristic length scales k of the fluid packets of both com- 
ponents are much longer than the domain of integratiorfl, if we 
just substitute /a, /b — * 5v, Av. Thus, from this analogy it is clear 
that we may set £ = / ve i as long as the Sobolev-like require- 
ment 6v » Cv t is satisfied, where Cv t is the velocity extent over 
which a photon of frequency x may be absorbed (that is, the 
velocity extent of a resonance zone). This limiting situation cor- 
responds to the case that the line profile can be represented by 
a delta function, so that the sharp edges of the resonance zones 
prevent any absorption at frequencies not Doppler shifted to the 
very line center, resulting in a localized radiative transfer. The 
optical depths in Eq. lA.ll are then understood to be the Sobolev 
ones. That is, r c] = r sm // ve i and r ic = T sm (x ic // vel ), with T sm the 
optical depth in the smooth case (Eq. [2j. 

However, especially in the outer wind (but, depending on the 
onset of clumping, also in the innermost wind, see Fig. |2]i we 
will generally have 6v < Cv t and the effective fraction of the ve- 
locity field over which photons can be absorbed by clumps will 
increase. The exact form of the radiation transport is then likely 
to be very complex. Nonetheless, let us in a first attempt try to 



4 We mention in passing that the Levermore et al. model also yields 
the result e~ s/ ' ! , with s the path length, for the normalized intensity in 
the limit of infinitely dense absorbers in a background vacuum. This 
is equivalent to the result for a fully porous wind obtained by, e.g., 
IOwockietal] ( l2004h . 



£ ~ /vei + Cr], (A.6) 

where rj = v t /Av is the effective escape ratio. Note the difference 
between this definition of r\ and that given in Paper! The two 
are related as 

7/=(l-/vel)A7old, (A.7) 

where denotes our earlier definition. The advantage of re- 
defining j] is that we may now separate out the porosity depen- 
dence in writing 

v t v t /(dvfj/dr) L r 

"=aT^— = T< (A ' 8) 

with h = 6r/f v - Ar - vpbt the porosity length of the medium 
in our geometry and L r = v t /(dvp/dr) the radial Sobolev length 
for a smooth velocity field. The coupling between vorosity and 
porosity becomes clear via 77. 

As defined, rj may in principle take arbitrarily high values, 
so for the examples in this paper we simply set £ = 1 when- 
ever £ > 1, because in a wind with a smooth velocity field the 
clumps obviously cannot absorb photons over a velocity space 
larger than that covered by the B velocity law. On the other hand, 
if we allow for clumps to be randomly positioned in velocity 
space, overlapping velocity spans will lead to a change in the 
effective coverage fractions. If velocity perturbations are suffi- 
ciently large, one may simply substitute £ — > (1 -e^) and permit 
£ to take arbitrarily high values. However, it is clear neither if ve- 
locity perturbations will be sufficiently large nor how to handle 
the case when more than one clump is crossed within a reso- 
nance zone. Thus we for now consider only the simple case of a 
smooth velocity field, deferring to future work a careful study of 
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these randomization effects. In any case, we note that our formal- 
ism recovers the smooth optical depth r sm in the limit Av « Cv t 
(as expected because then the individual clumps obviously are 
optically thin). 

Finally, Eq. lA.ll has the proper behavior in the limiting cases 
of a smooth or microclumped wind. For the former (x; c =l and 

R a = e- T "\ (A.9) 

and for the latter (r c i << 1) some simple algebra yields, 

R,*l-T sm , (A. 10) 

where we recall that this last result is expected because res- 
onance line formation depends linearly on the density (see 
Sect. ED, 



Recombination lines. The absorption part of recombination 
lines such as H„ may also be approximated as described above. 
Furthermore, since the source function in these lines can be pre- 
scribed (see Sect. 0, we can make a similar approximation for 
the re-emission part 

R em = S£(l - e^) +5(1 -£)(1 - e~ Ti <), (A. 11) 

where S is the source function at the resonance point in units of 
the continuum intensity. The total line profile R x is then given by 

= ^a,x + ^em,x- It is important to realize that the re-emission 
profile is much more influenced by non-radial photons than is the 
absorption part profile. Thus we replace the radial approximation 
for £ with a corresponding line-of-sight expression, £ — > ^ z , by 
substituting L, — > L z and h — » h/fj., where curvature effects for a 
clump have been neglected, and now obtain the final line profiles 
by performing the standard integration over a pre-specified num- 
ber of P-rays. The optical depths from the previous paragraph 
must be replaced by corresponding ones for recombination lines 
(see Sect. 13.21) . where of course care must be taken for the now 
angular dependent r. 

Also Eq. IA.1 II has the proper behavior for smooth as well 
as for microclumped winds. In the same manner as for the reso- 
nance lines, we obtain for the former 

R em = S(l-e- T ™), (A.12) 
and for the latter 

Rcm~ST sm f eh (A. 13) 

which is expected because recombination line formation de- 
pends on the square of the density (see Sect. 12. II ). 

Comparisons between the analytic approximations outlined 
here and numerical simulations using our stochastic wind mod- 
els and detailed radiative transfer codes are given in the main 
paper. 



